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ABSTRACT 

We present Roche tomograms of the secondary star in the dwarf nova system RU 
Pegasi derived from blue and red arm ISIS data taken on the 4.2-m Wilham Herschel 
Telescope. We have applied the entropy landscape technique to determine the system 
parameters and obtained component masses of Mi = 1.06 Mq, M2 = 0.96 M©, an 
orbital inclination angle of i = 43°, and an optimal systemic velocity of 7 = 7kms~^. 
These are in good agreement with previously published values. 

Our Roche tomograms of the secondary star show prominent irradiation of the 
inner Lagrangian point due to illumination by the disc and/or bright spot, which may 
have been enhanced as RU Peg was in outburst at the time of our observations. We find 
that this irradiation pattern is axi-symmetric and confined to regions of the star which 
have a direct view of the accretion regions. This is in contrast to previous attempts 
to map RU Peg which suggested that the irradiation pattern was non-symmetric and 
extended beyond the terminator. 

We also detect additional inhomogeneities in the surface distribution of stellar 
atomic absorption that we ascribe to the presence of a large star-spot. This spot 
is centred at a latitude of '^82° and covers approximately 4 per cent of the total 
surface area of the secondary. In keeping with the high latitude spots mapped on the 
cataclysmic variables AE Aqr and BV Cen, the spot on RU Peg also appears slightly 
shifted towards the trailing hemisphere of the star. Finally, we speculate that early 
mapping attempts which indicated non-symmetric irradiation patterns which extended 
beyond the terminator of CV donors could possibly be explained by a superposition 
of symmetric heating and a large spot. 

Key words: stars: novae, cataclysmic variables - stars: starspots - stars: late-type 
- stars: imaging - stars: individual: RU Peg - techniques: spectroscopic 



1 INTRODUCTION 



Cataclysmic variables (CVs - see lWarnerl (|l995h for an excel- 
lent review) are a class of interacting binary characterised by 
a white dwarf star, known as the primary, accreting material 
via an accretion disc or column from the secondary - typ- 
ically a late-type main-sequence dwarf star. CVs are often 
studied in order to gain insight into the processes of astro- 
physical accretion, but it is the secondary star that holds the 
key to unlocking our understanding of the creation and evo- 
lution of these systems and their corresponding behaviour. 
In the canonical scenario for the evolution of CVs, for ex- 
ample, the magnetic field of the secondary star is believed 
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to play a crucial role in draining away angular momentum 
from the system via a process known as magnetic braking 
(see, e .g., [Kraftl l|l967f ). iMestell (|l968l ) and ISpruit fc Rittej 
(1983^). This is thought to be responsible for maintaining 
mass transfer and for causing long period CVs to evolve to 
shorter orbital periods. At orbital periods of about 3 hours 
and stellar masses of around 0.25 M©, the secondary star 
is thought to become fully convective. Under the canoni- 
cal theory of CV evolution, this transition to a fully con- 
vective interior switches off the stellar dynamo and conse- 
quently mass transfer is quenched. Contact with the inner 
Lagrangian point is re-established at shorter periods through 
removal of angular momentum via gravitational radiation. 
This explains the so-called period gap in the period distribu- 
tion of the total CV population. Despite the acceptance of 
this model, there is little or no direct observational evidence 
for such an abrupt change in the angular momentum loss 
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mechanism and it is still unclear whether magnetic brak- 
ing is really the dominant loss mechanism above the period 
gap. Magnetic activity cycles have also been cited as possi- 
ble causes for a number of other observed CV phenomena, 
such as variations in orbital periods, lum inosity, mean out - 
burst durations and o utbu rst shapes fe.g. | Applegatel (|l992l ) , 
iRichman et all (|l994l ') and lAk et ai] (|200ll ')V Given the im- 
portance of understanding the underlying magnetic proper- 
ties of these stars, it is therefore highly desirable to try to 
obtain observational evidence for such fields, in the form of 
the sizes, distributions and lifetimes of star-spots. 

In addition, since CV secondaries are among the fastest 
rotating set of stars known, we may begin to use them as 
natural laboratories to provide tests of stellar dynamo the- 
ory, since they should display strong dynamo action. We 
may also begin to ask what effect tidal forces have on the 
emergence of magnetic flux tubes. It has been s uggested that 
tidal forces may suppress differential rotation (|Scharlemannl 
[1982) a nd force st ar-spots to form at preferred longitudes 
(Holzwa rth fc Schii ssler 2003). Since CV secondaries are 
heavily tidally distorted, it seems clear that these systems 
are perfect for probing these effects. 

However, despite the clear importance of magnetic be- 
haviour in the evolution of CVs and the potential gains to be 
had in our understanding of the so-called solar-stellar con- 
nection, unti l recently the r e has been little work undertaken 
in this field. IWebb et al] (|2002l ') used TiO bands to deter- 
mine a spot coverage factor of 22 per cent for the secondary 
in SS Cyg, although their work did not actually map the 
surface. Recently, however, the technique of Roche tomog- 
raphy has begun to allow the mapping of the surface of these 
stars, and results presented in papers III and IV of this series 
have reported the first detection of st ar-spots on a CV sec- 
ondary (AE Aqr - IWatson et al. 20061) followed by a similar 



mapping campaign for BV Cen I Watson et al. 2007af ). Both 
maps showed spots at almost all latitudes, as well as large 
near-polar spots similar to those found in Doppler imaging 
studies of single stars. In this work, we report on attempts 
to apply this technique to the dwarf nova RU Peg for the 
first time. RU Peg is a relatively bright (mv=13.1 rising to 
9.5 during outburst) U Geminorum type dwarf nova with an 
orbital period of 9 hours. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations 

Observations of RU Pegasi were made with the 4.2-m 
William Herschel Telescope (WHT) at the Roque de los 
Muchachos Observatory on the island of La Palma over the 
nights of 20-23 July 2004. A journal of these observations 
can be found in Table [T] The ISIS double beam spectro- 
graph, which makes use of a dichroic mirror to split the 
incoming light into blue and red channels, was used to ob- 
tain simultaneous blue and red spectra. The central wave- 
lengths for the blue and red spectra were selected to be 
4950A and 6400A respectively, yielding unvignetted ranges 
of 4550-5350A in the blue and 6070-6700A in the red. Both 
arms of ISIS made use of the R1200 gratings giving a disper- 
sion of about 0.23A/pixel. A two-pixel resolution of about 
28kms~^ in the blue and 22 km s~^ in the red was obtained 



for a slit width of 1 arcsec. Exposure times were selected 
in such a way as to balance the requirement for maximum 
signal while remaining short enough to avoid the effects of 
velocity smearing. Observations of the flux standard star 
SP 2148-1-248 were taken at or close to the zenith, in order 
to correct for the instrumental response and to remove any 
telluric absorption features in the red data. 

Since we require good relative flux calibration, we 
plac ed the nearby comparison star RU Peg Ref 2 
(Har rison et al.l 120041 ') on the slit to allow for slit-loss cor- 
rections and any variations in sky transparency. In addition 
to a number of observations of K and M spectral-type tem- 
plates, arc spectra were taken before and after each exposure 
(or multiple run of six ten minute exposures in the case of 
RU Peg) to enable the spectra to be wavelength calibrated 
and to correct for any flexing in the instrument as the tele- 
scope tracks the target. A set of tungsten flat fields and bias 
frames was also taken at the beginning and end of each night. 
The weather during the observing run varied from clear to 
heavy rain and fog. In effect, this resulted in the loss of one 
night out of four. Mechanical breakdowns also restricted the 
amount of observing we were able to accomplish. Amateur V 
band photometric observations (obtained by Roger Pickard 
and Gary Poyner of the British Astronomical Association 
Variable Star Section (BAA VSS), whom we gratefully ac- 
knowledge) showed that the system was not in quiescence 
during our observations. Rather, it was on the rise to peak 
outburst luminosity. The average blue and red arm spectra 
can be seen in Fig. [T] 



2.2 Data reduction 

The raw images were bias-subtracted using a median bias 
value obtained from the overscan region. Pixel-to-pixel vari- 
ations were removed by creating a master balance frame, 
generated from a set of tungsten flats, and dividing each 
raw frame through by the master balance frame. The spectra 
were finally sky subtracted and optimally extracted (jHornel 
[l98y) using the software package PAMELA. The Copper Ar- 
gon/Copper Neon (CuAr-l-CuNe) arc spectra were also op- 
timally extracted, from the same region of the chip that 
the object spectrum occupied. These were fitted with a 
fourth order polynomial giving an rms scatter of better than 
O.OOlA. The object spectra were then placed on this wave- 
length scale. The data were then flux calibrated using spec- 
tra of the featureless white dwarf star SP 2148-f 286. 

Roche tomography also requires that we correct for slit- 
losses and variable sky transmission. This is due to the 
rapidly varying and unknown spectral contribution of the 
accretion disc. Because ISIS is not a cross-dispersed spectro- 
graph, and so is capable of long-slit spectroscopy, we were 
a ble to place a non-va rying comparison star (RU Peg Ref 
2; [Harrison et al.ll2004h on the slit and take wide and nar- 
row slit exposures. We then multiplied each corresponding 
CV spectrum by the ratio of the wide slit comparison spec- 
trum to the narrow slit comparison spectrum. We should 
note that since the slit was not set to follow the parallactic 
angle, this is at the expense of some small effects arising 
from differential refraction. Since we treat both the red and 
blue arms independently during the Roche tomography re- 
constructions later, this reduces our exposure to such effects. 
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Figure 1. Left — WHT+ISIS average blue arm spectrum of RU Pegasi. Right - WHT+ISIS average red arm spectrum of RU Pegasi. In 
both cases the orbital motion of the secondary star has been removed prior to co-adding, and the continuum has been subtracted. 



Object 


Date 




UT Start 


UT End 


Tcxp (s) 


No. of spectra 


Comments 


BD+03 3568 


2004 


July- 


20 


23:08 


23:23 


300 


2 


G8V template 


BD+05 3625 








23:25 


23:35 


600 


1 


KOV template 


GAT 746 








23:52 


00:02 


600 


1 


KIV template 


RUPeg 








00:22 


02:11 


600 


13 


CV target 


SP 2148+286 








02:42 


02:52 


600 


1 


Spectrophotometric standard 


SP 2148+286 








02:53 


03:03 


600 


1 


Spectrophotomctric standard wide slit 


RUPeg 








03:15 


05:32 


600 


13 


CV target 


HD 154363 


2004 


July- 


21 


21:53 


21:55 


120 


1 


K5V template 


HD 181196 








22:03 


22:05 


180 


1 


K5V template 


HD 332601 








22:09 


22:14 


600 


1 


K5V template 


HD 160964 








22:24 


22:28 


240 


1 


K4V template 


SP 2148+286 








00:34 


00:44 


600 


1 


Spectrophotometric standard 


SP 2148+286 








00:45 


01:06 


600 


2 


Spectrophotomctric standard wide slit 


RUPeg 








01:28 


05:41 


600 


45 


CV target 


BD+00 3459 


2004 


July- 


22 


22:49 


22:50 


60 


1 


K3V template 


RUPeg 








02:35 


04:37 


600 


22 


CV target 


BD+30 3627 








04:47 


04:52 


300 


1 


K2V template 


SP 2148+286 








04:56 


05:06 


600 


1 


Spectrophotometric standard 


HD 333388 








05:14 


05:16 


150 


1 


K8V template 


BD+01 3236 


2004 


July- 


23 


21:39 


21:47 


450 


1 


K7V template 


GJ 9606 








22:10 


22:20 


600 


1 


MOV template 


GJ 9645 








22:27 


22:35 


450 


1 


K6V template 


SP 2148+286 








22:41 


23:15 


600 


3 


Spectrophotometric standard 


RUPeg 








02:26 


03:48 


600 


13 


CV target 



Table 1. Journal of observations for the July 2004 WHT/ISIS run. The first column gives the object name, followed by the date and 
start/end times of the observations. The final three columns indicate the exposure times (for the blue arm of ISIS), the number of spectra 
obtained and an indication of the type of spectra taken. 



3 THE SECONDARY STAR 
3.1 Projected rotational velocity 

In order for accurate orbital phases to be calculated for our 
spectra, a new orbital ephemeris was required (see section 
13. 2|) . However, in order to provide the most accurate cross- 
correlation between our target spectra and the template, we 
also need an estimate for the mean line width of the tar- 
get spectrum. To be able to measure the mean projected 
equatorial rotational velocity of the secondary, t;sini, we 
require an average spectrum for maximum signal to noise. 
Simply creating a mean spectrum at this point would be fu- 



tile, as the orbital motion of the secondary would effectively 
smear out the absorption line profiles due to the Doppler 
shifts throughout the orbit. We therefore subtracted the ra- 
dial velocity of each spectrum individually, having deter- 
mined their velocities through cross-correlation with a high 
signal-to-noise spectral-type template star. The spectra are, 
of course, composites of all the components in the CV. To en- 
sure the secondary star line strengths were preserved across 
the region of interest, the spectra were first rebinned onto 
the same velocity scale using sine-function interpolation and 
then normalised by dividing through by a constant. Finally, 
the spectral continua were fitted with a third-order spline 
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Figure 2. x^ versus i;sini obtained through cross correlation 
with a K5V template. The minimum is located at a projected 
rotational velocity of 88kms~^. Note (see text) that the broad- 
ening increment of 5 km s~^ has been reduced to 1 km s~^ around 
the minimum to allow its position to be better determined. 



which was then subtracted. An identical procedure was ap- 
plied to the spectral type template with the added compli- 
cation that the spectral type standard was corrected for its 
systemic velocity. 

The majority of our template stars do not have prior 
systemic velocity measurements and so we used Least 
Squares Deconvolution (LSD - see section |4]) to calculate 
a mean absorption line profile for each spectrum. This mean 
profile was then fitted with a gaussian and the location of the 
peak of the fit was taken as the template's radial velocity. 
Table 2 lists the measured radial velocities of our spectral 
type templates. A previous estimate of the r adial velocity 
of our K5V template GL 653 can be found in lEvand (119671 ) 
who reported a value of 33.6 ±2 km s~^. This is within 1.5cr 
of our LSD derived estimate of 30.6 km s^^. 

A first iteration on the radial velocity curve was ob- 
tained by cross-correlating the RU Peg spectra with a K5V 
template that had been rotationally broadened by an arbi- 
trary amount, lOOkms"^ in this instance. We used the red 
arm data, specifically wavelengths between 6200A & 6500A, 
since this range contains a number of relatively strong stel- 
lar absorption lines and also some kno wn temperature and 
gravi ty sensitive lines for K type stars (IStrassmeier fc Fekell 
I1990I ). This first pass made it possible to subtract approxi- 
mate radial velocities of the spectra and create an average 
RU Peg spectrum. 

The template spectrum was then artificially broadened 
by 5kms~^ increments, multiplied by some constant and 
then optimally subtracted from the average RU Peg spec- 
trum. In this process, the fraction of the template that best 
removes the secondary star lines is calculated using an opti- 
misation technique that minimises the scatter between the 
template and the target in a process known as optima l sub- 
traction (a term first used by iDhillon fc MarshI l|l993l ) and 
widely used since in the context of secondary stars in CVs). 
Fig. |2] shows the results of this method for the K5V tem- 
plate, where we have decreased the broadening increment to 
lkms~^ around the curve's minimum. The optimum value 
of vsini is 88kms^^. 



K type subclass 



Figure 3. RU Peg secondary star spectral type. The x-axis refers 
to the subclass of K dwarf template used in the optimal subtrac- 
tion procedure described in section [XT] 



To test whether our choice of spectral type template in- 
troduces systematic errors, or indeed if there is any trend in 
the measured value of v sin i with spectral type, the proce- 
dure has been repeated for all our spectral type templates. 
We find that the choice of spectral type does not significantly 
affect the resulting optimal value for v sin i - the standard 
error on the mean of our sample is only 0.5kms~^. In each 
application of artificial broadening, we have assumed a sim- 
ple linear limb darkening law with a coefficient of 0.5, 



7(0) 



= l-a(l-/i) 



(1) 



where /i = cos 5 {5 being the angle between the line of sight 
and the emerging intensity), 1(0) is the specific intensity 
at the centre of the stell ar disc and a is the linear limb 
darkening coefficient (see IClareJ 1 19981 ). Again, to test for 
the possibility of introducing systematic errors, the entire 
procedure has been run using a range of values from to 

I for the limb darkening coefficient. The optimum values of 

II sin i did not differ from each other by more than 3kms^^. 
We therefore adopt vsini = 88 ± 3kms~^. 

The spectral type of the secondary star can also be esti- 
mated using optimal subtraction, given a full set of spectral 
type templates. Fig.Oshows the reduced x^ against K dwarf 
subtype which displays a mininmm at K5V. Again, we tested 
the effect limb darkening has on the outcome. We found no 
change in the optimum spectral type while varying the lin- 
ear limb darkening coefficient. Our estimated spectral- type 
compares well wi th the previous esti mated spectral type of 
K3V a dopte d by iFriend et all (|l990l ) following the work of 
IWadd (|1982| ). 



3.2 Radial velocity curve and ephemeris 

The most recent orbital ephemeris for RU Peg was de- 
rived from second ary star radial velocity variations by 
iFriend et al.l l|l990l ). Given the long gap between their obser- 
vations and the ones reported here, a large error may have 
accumulated in the orbital zero point. We have derived a 
new ephemeris by cross-correlating our red arm RU Peg data 
with the best-fitting K5V spectral type template GL 653 to 
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produce a series of cross-correlation functions (CCFs). The 
RU Peg and the GL 653 spectra are treated in the same 
way as described in section [371] with the template being ar- 
tificially broadened by the optimum amount found through 
optimal subtraction, in this case 88kms~^. In addition, the 
radial velocity of the K5V template has been removed. A 
cross correlation function was generated for each RU Peg 
spectrum, and the radial velocity curve obtained by fitting 
a sinusoid through the CCF peaks. The following ephemeris: 



To = 2453207.53376 ± 0.00009, 



(2) 



was obtained where the orbital period has been fixed at 
0.3746 days (|Stoveij[T98il '). and To is the time of inferior 
conjunction of the secondary star, in the heliocentric ref- 
erence frame. The quoted uncertainty is statistical, with 
no allowance for any possible systematic effects. This new 
zero point has been used in all subsequent analysis. The 
robustness of this method of generating a radial velocity 
curve through cross-correlation against an incorrect tem- 
plate choice has been tested. A series of CCFs was produced 
using the full range of K dwarf templates at our disposal, and 
it was found that, providing the templates were correctly 
broadened, there was no significant change in the value of 
the semi-amplitude or systemic velocity, and therefore the 
zero- crossing. 

In Fig. |4l the radial velocity points have been assigned 
an orbital phase based on the new ephemeris, and fitted with 
a circular fit of the form. 



V {(j)) = 'Y + K2 sin 



2tt {t - to) 



(3) 



where K2 is the semi-amplitude of the secondary star, 7 is 
the systemic velocity, V {(j)) is the radial velocity at a phase 
<j), t is the time of observation, to is the zero crossing and 
Porb is the orbital period. From this fit we obtain a value 
for the radial velocity semi-amplitude of 128.2 ± 0.6kms~^. 
This com p ares with a value of 111 ± 8kms~^ fou nd in 
iKiplingeil (|l979l ). 121 ± 2 kms~^ from IStoveJ (|l98ll ) and 
121 ± 2kms"^ derived in I Friend et~all (|l990l '). We should 
add at this point that there is a clear systematic trend 
in the residuals of the circular fit of Fig. 3] This is most 
likely the result of irradiation, making such conventional 
studies susceptible to systematic errors since the centre-of- 
light of the secondary no longer represents its centre-of-mass. 
Our value for the s ystemic ve l ocity of 10 ± 0.5kms~^ com- 
pares well with [Friend et al.l ([1990|) who report a value of 
13 ± 2 kms~^, and Stov er's value of 5 ± 4kms~^. The anal- 
ysis of iKiplingeiJ (| 19791 ) yielded a slightly more discrepant 
va lue of — l±5k ms~^ using the radial velocity data obtained 
in iKrafj ([1962j, though the value changes to 18 ± 6kms~^ 
when measuring 7 from a disc emission line radial velocity 
curve. 

Additional motivation for measuring the projected ro- 
tational velocity and K2 is that one can determine a value 
for the mass ratio and in turn, the se mi-amplitude of the 
primary star. From lHorne et al.l (|l986r ). the predicted rota- 
tional velocity for a Roche-lobe filling star is given by: 



i?L 



K2 (1 -f g) . 



(4) 




Comb ining this with the approximation of lSmith fc DhillorJ 
l|l998h for the volume- averaged Roche lobe radius (assumed 



Figure 4. RU Peg (ISIS red arm) radial velocity data (filled 
circles) and circular orbital fit. 



equal to the radius of the secondary star) gives the following 
relation: 

vs\ui = QA7K2q^''-^{l + qf'^ . (5) 

Rearranging for g, the mass ratio (defined as M2/M1), and 
inserting our values for vsini and K2 derived earlier, equa- 
tion ([5} can be solved to obtain g = 0.88rb 0.02. This is h igher 
than the values of 0.78 and 0.79 found in lStoverl (|l98ll ) and 
iKiplingeiJ (| 19791 ) . both obtained directly through measure- 
ments of both Ki and K2 ■ Given the relatively poor quality 
of their emission line radial velocity curves, it may be argued 
that the value derived here from accurate rotational velocity 
measurements is somewhat more reliable, although our value 
for K2 has not fully ta ken into account the residuals in the 
fit. lFriend et al.1 (|l990l ) derive a value of g = 0.82±0.15 from 
their rotational broadening measurements that is slightly 
closer to ours. However, the value presented here is con- 
sidered an improvement given the higher resolution of our 
data. Our mass ratio and K2 imply a primary star radial ve- 
locity semi-amplitude Ki = 112.6 ± 0.25 kms"'^. iKiplingerl 
(| 19791) reports a value of 88 ± lOkms^^, significantly lower 
than our value, from direct measurement of the Ha emis- 
sion wings; althoug h the fi t is rea sonable, the reported rms 
error is 26kms~^. IStoven (|l98ll ). using a similar method 
for measuring the emission wing velocities, finds a value of 
94 ± 3kms~'^ fitting a circular orbit to the radial velocity 
data (with an rms error of 12kms~^). We note here that no 
attempt has been made to directly measure the emission line 
radial velocities in this work due to the complicated nature 
of the emission line profiles. Since the system is in outburst, 
we see that H/3 and Ha have two anti-phased components - 
one from the disc and one from the irradiated surface of the 
secondary. 



4 LEAST SQUARES DECONVOLUTION 

While single absorption-line studies of CV donor stars 
using intermediate resolution spectrographs on 4-m class 
telescopes have been successful in mapping large-scale ir- 
radiation patterns on their surfaces (e.g. IWatson et al.l 
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Spectral Type 


Name 


Radial Velocity ( km s ^) 


G8V 


BD+03 3568 


-0.3 


KOV 


BD+05 3625 


-31.9 


KIV 


GAT 746 


-50.0 


K2V 


BD+30 3627 


-57.3 


K3V 


BD+00 3459 


-5.3 


K4V 


HD 160964 


-20.8 


K5V 


GL 653 


30.6 


K8V 


HD 333388 


-7.1 



Table 2. Heliocentric radial velocities of spectral type templates 
measured via a Gaussian fit to the mean profile. 



[2003), they have not been successful at mapping star- 
spots because of th e limited res ol ution and signal-to- 
noise of such data. IWatson et al.l l|2006l ). however, ap- 
plied the technique of Least Squ ares Deconvolution (LSD - 
iDonati k. Collier Cameronl ll997l to produce a mean profile 
by combining several thousand absorption lines from the sec- 
ondary star in AE Aqr observed with the UES echelle spec- 
trograph on the WHT. This achieved the required signal-to- 
noise to observe the effects of star-spots on the line profile. 

LSD requires the positions and strengths of the atomic 
absorption lines present in the stellar spectrum to be known. 
In this work, we use an atomic line li st generated by the 
Vien na Atomic Line Database (VALD. iKupka et al.l l|l999l . 
120001 )). VALD offers a simple online interface for the selec- 
tion of the required data; all that is required are values for 
the stellar effective temperature and surface gravity, the rele- 
vant wavelength range, micro-turbulence (due to small scale 
motions of convective cells at the stellar surface, taken here 
to be «lkms~^) and a detection limit for the strengths of 
the normalised lines. The detection limit in Doppler imag- 
ing applications is often taken to be about 0.4 (this is the 
lower limit to the central line depth, ignoring line blending 
and expressed as a fraction of the continuum) . Since we are 
dealing with intermediate resolution data rather than echelle 
data, we have used a more conservative 0.5 (this may still 
seem optimistic; however tests described below indicate that 
the input line list has little effect on the shape of the mean 
profile). The spectral type of RU Peg has been estimated in 
this work as K5V - the closest model available in the VALD 
database corresponds to a stellar atmosphere with an effec- 
tive temperature, Toff, of 4600K and log {g/ cms~^) — 4.6. 

The composite nature of the CV spectrum complicates 
the situation with regard to treating the spectral continuum. 
There is a varying and unknown contribution from the ac- 
cretion regions, as well as a changing continuum slope due to 
any time variability and the changing sky projection of the 
accretion disc and hot spot. Th is means that the creation 
of a master continuum fit (e.g. ICollier-Cameron fc Unruhl 
[l99J) is impossible and we are forced to remove the contin- 
uum on a spectrum by spectrum basis. Simply normalising 
the spectra is also inappropriate, as the varying secondary 
star contribution would result in the relative line strength 
changing from one spectrum to the next. We therefore con- 
tinuum subtracted the data (after flux calibration and slit- 
loss correction) by subtracting a smoothly varying fit for 
each spectrum. This fit is created by placing spline knots at 
positions that are relatively free of absorption lines, while 
simultaneously ignoring strong disc emission lines. 



This method of fitting the continuum may appear rather 
subjective, since it is often difficult to decide exactly how the 
continuum varies across a spectrum. Tests were carried out 
to determine the effects of incorrect continuum fits on the 
final mean profile by constructing two continuum fits to the 
same spectrum - the first being our best guess at the true 
continuum and the second was a fit which was clearly wrong. 
We find that the continuum fit does not have a significant 
effect on the shape of the resultant mean profile. However the 
continuum fit chosen does, of course, have an impact on the 
level of the continuum seen in the mean profiles. In general, 
it was found that the final continuum fit selected was at 
too high a level, resulting in the (subtracted) continuum 
wings of the mean profiles lying well below zero. This is 
straightforward to remedy by iteratively lowering the level 
of the continuum fit until the desired zero continuum level 
is achieved in the LSD profile. Again, this process was found 
not to alter the shape of the mean profile. 

It should be noted that since our line list contains nor- 
malised line depths and our spectroscopic data are contin- 
uum subtracted, a further step was needed to prepare the 
line list for use with our data. Each line depth in the line 
list was multiplied by the value of a continuum fit to a K5V 
spectral type template at that wavelength. In this way, the 
line list now contains a set of lines that have the correct 
relative depths and may be used with continuum subtracted 
data. The sensitivity of the LSD process to the choice of line 
list used has been investigated by carrying out LSD on the 
same RU Peg spectrum using a range of spectral type line 
lists t hat cov e rs the spectral types KOV to K8V inclusive. 
As in iBarned (|l999l ). there is no evidence that slight errors 
in the assumed spectral type, and therefore an incorrect line 
list, introduces any appreciable systematic errors in the final 
average profile. 



5 ROCHE TOMOGRAPHY 

Roche tomography is an astro-tomographic technique de- 
veloped to image t he surface of the secondary star in cat - 
aclysmic variables JRutten fc DhillonI [l993 : 1 Watson! 120021 '). 
The procedure is, by now, robustly tested and arte- 
facts caused__b^;_s^;2teinatic errors are well understood 
fsee IWatson fc Dhflloiil I2OOII for details). Early work fo- 
cused on map ping single absorption or emission lines (see 
I Watson et al.l 12003) while, more recent ly, efforts have been 
made to map star-spot patterns (see IWatson et al.l |2006| . 
2007a). 

Roche tomography takes as its input a trail of phase 
resolved spectral line profiles, and produces a map of the 
secondary star in three-dimensional real space. The key ad- 
vantage this technique has over other indirect imaging meth- 
ods (such as Doppler tomography) is that it makes use of all 
the information contained in the line (i.e. its fiux, width and 
Doppler shift) to map the star. Roche tomography is, how- 
ever, very similar to Doppler imaging of single stars, but for 
the exceptions that in Roche tomography the star is Roche- 
lobe filling, rotates around the binary centre-of-mass, and 
is tidally locked. Roche tomography also ignores the contin- 
uum. Successfully applying the technique is often somewhat 
harder than Doppler imaging. Roche tomography targets are 
generally much fainter than the targets of choice in single 
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Figure 5. Roche input trailed spectrum, computed data and residuals for RU Peg blue data. Dark grey scales depict deeper line depths 
relative to lighter grey scales. 



star mapping programs, resulting in spectral profiles of lower 
signal-to-noise and poorer quality maps. 



We use t he tomography code Rochey developed by 
_Watson (^200 j) and f irst applied with the addition of LSD by 
Watson et alj (|2006l ). The synthetic LSD profiles were con- 
structed by assuming that each tile contributes a Gaussian 
profile to each line, scaled by projected area and intensity. 
Although temperature and abundance will affect the profile, 
so that the spot contribution to the profile will be slightly 
different from one line to another, we don't have enough 
information to take that into account and we simply com- 
pute an average profile for each tile, and use that average 
profile in the mapping. The LSD process is not very sen - 
sitive to temperature and metallicity effects (|Barneg|l999h . 
Although for single stars it is possible to assume for each tile 
some combination of 'spot' and 'immaculate' assumed spec- 
trum, that is not possible here because irradiation means 
that there is no single 'immaculate' spectrum - and indeed 
the contribution of irradiation is something that we want to 
obtain. We recognise that this means that what the Rochey 
code produces is essentially a contrast map. 



5.1 System parameters 

The use of incorrect system parameters in the Roche tomog- 
raphy reconstru ction results in artefacts being mapped onto 
the stellar grid (| Watson fc Dhillonll2001h . For example, a set 
of incorrect values for the masses and the orbital inclination 
will result in the code not being able to properly fit the varia- 
tions of the profile widths along the trail. Roche tomography 
allows a natural method of selecting optimal values for these 
parameters, since only the correct values will give the map 
of maximum entropy. We carry out reconstructions for a set 
of component mass pairs, fitting to the same x^ each time, 
until a full grid of values is obtaine d - the so-called entropy 
landscape (|Rutten fc DhillonllTooi '). 



Limb Darkening 

An incorrect treatment of the limb darkening for the sec- 
ondary star may also result in extra structure being mapped 
onto its surface. In this work we make use of a square-root 
limb darkening law of the following form (|Claretlll998l ). 



/ ifi) = Jo [1 - aA (1 - m) - bx (1 - Vm)] 



(6) 



where ax and ^a are the limb darkening coefficients at a 
wavelength A. IClaretl (119981 ) emphasises the non-linearity of 
limb darkening laws applied to model atmospheres of late K 
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and M type main sequence stars; for the coolest stars, the 
departure from hnearity is drastic. 

Since hmb darkening is a function of wavelength, we de- 
termine an effective central wavelength for the blue arm ISIS 
data of 4875A that takes into accoun t the line depths and 
spectral noise and is calculated after IWatson et alj (|2006l ) 



also true here, where for the blue data, the optimum sys- 
temic velocity of 7kms~^ is constant for the range of incli- 



nations probed. Fig. 6(c) shows a slice through the entropy 



usmg. 



'^ccn — 






(7) 



where di and cri are the flux and error at wavelength Ai, 
respectively. 

We assume a secondary star spectral type of K5V (see 
section I3.1|l and adopt values for the effective tempe rature 
of 45 57K and log surface gravity (in cms" ) of 4.6 (JGravl 
Il992r ). The closest model available in the dat abase of non - 
linear limb darkening coefficients provided bv lClaretl (|l998l ) 
has an effective temperature of 4660K and surface gravity of 
4.6. We obtained values for the limb darkening coefficients 
in the Johnson B and V bands and linearly interpolated 
across them to arrive at values corresponding to the cal- 
culated efl'ective central wavelength of our spectra. These 
coefficients, a\ — 0.8062 and bx = 0.0538, are used in all 
subsequent reconstructions. For the red arm ISIS data with 
an effective central wavelength of 6426A, we obtain coeffi- 
cients of ax = 0.482 and bx ~ 0.347 for the Johnson R band 
(Accn ~ 6300A). Note that it was not possible to interpolate 
values of the coefficients at our central red wavelength, be- 
cause of a lack of limb darken ing data for photometric bands 
later than R in lClaretl (|l998l '). 

Note that we have used continuum limb-darkening co- 
efficients. We recognise that coefficients for the cores of ab- 
sorption lines are considerably smaller; however, the pre- 
cise values for line limb-darkening coefficients depend on 
many factors, such as the angle of incidence, and there- 
fore require detailed and time-consuming calculations (e.g. 
ICollins II fc Truajg (| 19951 ')'). We have not attempted these, 
for three reasons. Firstly, it is the wings of the lines, where 
the difference from the continuum is least, that constrain 
the rotation and therefore the mass ratio. Secondly, we have 
tried varying the limb-darkening coefficients, and even with 
no limb darkening the maps still show the same gross fea- 
tures, although the parameters change slightly. Thirdly, the 
coefficients we have chosen give a somewhat lower x^ than 
other values we tried, even though the results are not very 
sensitive to the choice. 



Systemic velocity, inclination and component masses 

Figs 6(a) and (b) show surfaces of maximum entropy for 
a set of inclinations and systemic velocities for the blue 
and red ISIS data, respectively. For each systemic veloc- 
ity/inclination pair an entropy landscape has been con- 
structed, and the maximum entropy value extracted. In 
practice, an inclination value is first selected based on pre- 
vious estimates and a series of entropy landscapes are gen- 
erated for a range of systemic velocities. 

Previous Roche tomography work (Watson et al. 2003, 
2006, 2007a) has shown that for almost any choice of in- 
clination (and indeed for a wide range of mass pairs) the 
entropy landscape yields the same systemic velocity. This is 



surface shown in Fig. 6(a) at an inclination of 43°, show- 
ing a peak in the entropy at 7kms~^. Similarly, Fig. 6(d) 



shows a cut through the red maximum entropy surface at 
an inclination of 42°, where the entropy clearly peaks at a 
systemic velocity of 6.5kms~^. It is comforting to note that 
the optimal systemic velocity derived from two essentially 
independent data sets is constant for the full range of incli- 
nations probed, and has essentially the same value for the 
two data sets. 

We are unable to assign a fully rigorous (Monte Carlo) 
error determination as such a procedure would be vastly 
computationally expensive (see also below). Test reconstruc- 
tions, using values for the systemic velocity on either side of 
our optimum value, result in rapid degradation of the image 
quality, as also found by Watson et al. (2007a), who reported 
that "reconstructions were almost impossible for assumed 
systemic velocities that differed by more than ±2kms~^ 
from the optimal value". In addition, we found that viable 
reconstructions are impossible at inclination values ±10° 
from our optimum value. It is therefore unlikely that the 
errors exceed half these extreme values. 

While our optimal systemic velocity is close to the value 
derived from the more traditional radial velocity study in 
Section [31 which yielded 10.2 ± 0.4kms~^ for the systemic 
velocity, it does not agree with it within the error bars. This 
is not too surprising, since the value obtained through the 
radial velocity curve is subject to an incomplete treatment 
of the surface distribution of flux on the secondary, which 
is known to cause systematic errors in the systemic velocity 
determination. Indeed, as noted earlier, there remain resid- 
ual radial velocities of up to ~10% at some orbital phases 
afte r fltting a circula r orbit to the radial velocity data. 

iKiplingen (|l979l ) used the earlier radial velocity data of 
iKraftI ((1962!) to derive a value of — 1 ± 5 kms~^ for the sys- 
temic velocity (though he was unable to detect the secondary 
star in this work due to severe veiling of the secondary lines 
by a bright disc o n decline from outburst). It should be noted 
that the work of lKrafj l| 19621 ) represents the earliest detec- 
tion of RU Peg as a binary, and the radial velocity data re- 
ported there are derived from very low resolution data, with 
the subsequent or bital fit showi ng a large scatter. Our value 
is consistent with lStoveii (|l98l|), who reports a systemic ve- 
locity of 5±4 km s~^ derived from observations of the system 
taken in quiescence with higher quality data. iFriend et al.l 
(1990) quote a slightly larger value of 13 ± 2kms~^ using 
the Nal doublet, though they question whether or not the 
secondary star semi-amplitude is truly accurate, or if it is 
biased by low level irradiation. It is noted here that an ac- 
curate value for the radial velocity of the template used in 
the cross-correlation is also important, since this must be 
subtracted from the CV radial velocity data prior to orbital 
fitting. In our case, this has been accurately measured using 
LSD. 

Once the systemic velocity has been determined, it is 
then possible for the inclination and the component masses 
to be determined. These are solved for essentially simulta- 
neously, although for computational convenience we in prac- 
tice chose a value of i and produced an entropy landscape 
for the two masses, then repeated this process for a range 
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inclination systemic velocity (km/s) 

(a) Surface of maximum entropy (blue ISIS data). 



in'='i'^^^i°" systemic uelocity (km/s) 

(b) Surface of maximum entropy (red ISIS data) . 




(c) Maximum entropy as a function of systemic ve- 
locity (ISIS blue) 



systemic velocity {km/s} 

(d) Maximum entropy as a function of systemic ve- 
locity (ISIS red) 





(e) Maximum entropy as a function of inclination 
(ISIS blue) 



(f) Maximum entropy as a function of inclination 
(ISIS red) 



Figure 6. Optimal system parameters. In the top panels, the value of each point on the surface is obtained from an entropy landscape 
carried out at a particular combination of orbital inclination and systemic velocity. The optimum systemic velocity is clearly independent 
of the assumed inclination. On this scale however, it is unclear where the peak inclination lies. The lower panels show sections through 
the surfaces (a) and (b), and allow the peaks to be determined (see text). 



of i, choosing finally the inclination and mass pair with the 
highest overall entropy. 

We first show how the entropy depends on the choice 



of i. Fig. 6(e) shows the maximum entropy values extracted 
from a series of mass landscapes obtained over a range of 
inclinations with a fixed systemic velocity of 7kms~^ using 



our blue data. Similarly in Fig. 6(f) the maximum entropy is 
plotted assuming the systemic velocity of 6.5kms~^ derived 
earlier, for the red ISIS data. For the blue arm, we obtain a 
value of 43° For the red arm data, this value shifts slightly 



to 42°; just as we chose a unique value of 7 based upon the 
relative quality of the mean profiles obtained from the blue 
and red data, we also chose to adopt the blue arm derived 
value of 43° for the inclination of RU Peg. 

It is at first sight surprising that an inclination can be 
deduced for a non-eclipsing system. However, the inclination 
affects the line profile in two ways. First, the variation in the 
apparent width of the stellar absorption line-profile around 
the orbit depends on i (e.g. lShahbaz 199^ ); second, the way 
the spot bumps traverse the profile depends quite strongly 
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Data 



K2 (kms-l) Ki (kms-l) 



Red arm, from tomogram 


131.2 


116.1 


Blue arm, from tomogram 


133.5 


123.6 


Direct, red arm, this paper 


128.2 


112.6 


Friend et al. (1990) 


121±2 




Stover (1981) 




94±3 


Kiolineer fl979) 




88±10 



Table 3. A comparison of the radial velocity semi-amplitudes im- 
plied by our maximum entropy results with direct measurements. 



on i (for example, there would be no motion of the bump 
at all if the inclination were zero). The maximum entropy 
method disentangles this dependence and yields an inclina- 
tion. As with the systemic velocity, we would ideally like to 
estimate all uncertainties rigorously by undertaking a full 
Monte Carlo study; however, we have estimated that this 
might take up to a month of computer time, and no-one do- 
ing Roche tomography studies has yet attempted a full error 

ana lysis. 

Friend et al.1 l|l990l ) find a value of 33° ± 5° for the in- 



clination derived by assuming a main-sequence mass-radius 
relation for the secondary. However, if the primary sta r is a 
carbo n-oxygen white dwarf (i.e. Mi — 0.55-1.2 M© - llbeni 
1 19851 ). then the allowed range of inclinations for RU Peg is 
34°- 48°, which our value lies well within. The assumption 
that the secondary obeys a main sequence mass radius rela- 
tion is also applied in iKiplingeiJ (| 19791 ) who quotes a lower 
limit of 30° . Our valu e is, however, in good agreement with 
ISion fc UrbanI (|2002l ). who determine a value of 41° from 
model fitting of the disc spectrum. 

The mass entropy landscapes corresponding to the 
above values of systemic velocity and inclination are shown 
in Fig. [T] From these, we derive estimates for the component 
masses of Mi = 1.08 M© and M2 = 1.0 M© for the blue 
data, and Mi = 1.04 Mq and M2 = 0.92 M© for the red 
data. While the derived primary star mass values are quite 
similar, there is a greater difference between the values for 
the secondary star masses. We see no reason to preferentially 
adopt one estimate of the mass (either derived from the blue 
or the red data) over the other. Instead we adopt an average 
of these values and estimated uncertainties of ±0.04 M© for 
the primary masses and ±0.08 M© for the secondary masses 
which take into account the scatter between the values, but 
not any possible systematic effects; the real uncertainty is 
very likely to be larger. 

Having determined the best-fit system parameters, it 
is of interest to see what radial-velocity amplitudes are im- 
plied by these values and how they compare with previous 
direct determinations. Because the system parameters from 
the blue and red data differ, we consider these data sets 
separately. The results are given in Tabled where they are 
compared with our direct values from Section [3. 2 1 and values 
from the literature. 

The deduced K2 values from the two independent 
datasets agree satisfactorily and are not very different from 
the direct value found in Section [T2l More surprisingly, they 
are slighter larger than the direct value, which is not what 
would be expected if the difference were due to irradiation. 
This probably means that irradiation has no dramatic ef- 
fect on the radial velocity amplitude in this long-period. 



low inclination system, and that the values are consistent, 
implying that the uncertainties on each are about 4kms~^. 
While our K2 value of ~130kms~ ^ differs significantly from 
that found bv lFriend et al.l (|l990l ). we note that their data 
was of significantly lower resolution (at best 75kms~^) and 
this may account for the difference. The deduced Ki val- 
ues are rather more different from one another. If they are 
also assumed to be free from systematic errors, then this 
would indicate an uncertainty in the derived Ki values of 
~5kms~^. However, the three values from this paper are 
sufficiently si milar t o rule out t he conside r ably smaller val- 
ues found bv lStoved l|l98ll ) and lKiplingeil (|l97'9l ). 



5.2 Surface images 

Fig. [5] shows the input trailed spectrum of the blue-arm 
LSD absorption profiles together with the computed data 
from Roche tomography and the corresponding residuals. 
The Roche tomogram produced from this trail is shown in 
full in the left hand panel of Fig. [S] Perhaps the most visually 
striking feature is the large, dark, high latitude star-spot. 
This is a surprising result given the intermediate resolution 
of the ISIS spectrograph which we had assumed would be 
insufficient to map spots. However, given the inclination of 
the system (43°), polar features would be visible through- 
out the whole orbit and would contribute at every phase, 
thereby lending themselves to detection with even relatively 
low resolution instruments. 

Similar large high latitude spots are commonly found in 
Doppler imaging studies of rapidly rotating single stars and 
have also been recent ly detected in tw o other CVs - AE Aq r 
(|Watson et al.1 I2OO6I ) and BV Cen (| Watson et al. 2007al ). 
These features are in clear contrast to sun-spots, whose lo- 
cations are rarely more than 30° from the solar equator. 
In the case of RU Peg, we have determined that the cen- 
tre of the spot lies at a latitude of ~82°, and covers an area 
approximately 4 per cent of the total surface area of the sec- 
ondary star. Wh at causes the emergence of such high lat- 
itude star-spots? ISchiissler fc Solankil l| 19921 ) propose that 
the strong Coriolis forces dominate over the buoyancy of 
emerging magnetic flux tubes from deep within the convec- 
tive zone and force them to follow a path almost parallel 
to the rotational axis. Although this conclusion was drawn 
for single stars, the dynamics of magnetic flux emergence is 
essentially the same for synchronously rotating stars. An al- 
ter native_scenario_of migrating star-spots has been suggested 
bv lVogt et al.l (|l999l ). who report on the apparent clockwise 
spiral migration of small long-lived spots towards the pole, 
perhaps merging with the polar spot. Whether these mi- 
grating spots reinforce or cancel out the polar spot is of 
course dependent on their (as yet) unknown polarity. The 
star-spot in our map also appears to be deflected slightly 
towards the trailing edge of the star. This apparent deflec- 
tion is also clear in the Roche tomograms of AE Aqr and 
BV Cen by Watson et al. (2006, 2007a). A comparison of 
the Roche tomogram of RU Peg to those of AE Aqr and BV 
Cen is s hown in Fig.[5]along wi th a preliminary map of V426 
Oph bv I Watson et al. f2007b)l . Is this a common feature of 
tidally distorted rapidly rotating stars? If so, what is the ex- 
act mechanism that produces it? Perhaps it is caused by the 
orbital motion of the binary. Interestingly, there is no such 
high-latitude spot present in the V426 Oph map, though 
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Figure 7. Le/i - Entropy landscape for the RU Peg blue data, assuming an inclination of 43° and a systemic velocity of 7kms^^. 
Black pixels represent mass combinations for which no acceptable solution could be found. The black cross indicates the entropy peak, 
corresponding to primary and secondary masses of Mi = 1.08 Mq and M2 = 1.0 Mq. Right - Entropy landscape for the RU Peg red 
data, assuming an inclination of 42° and a systemic velocity of 6.5kms~^. The component masses giving rise to the maximum entropy 
are Mi = 1.O4M0 and M2 = 0.92 Mq. 



this system was caught in an extended (and unusual) qui- 
escent state, where its regular outbursting behaviour had 
essentially halted for over a year. Thus, if magnetic activ- 
ity is responsible for the behaviour of CVs then the map 
of V426 Oph may not be representative of the system as it 
would normally appear. 

Given the resolution of our RU Peg data, this is prob- 
ably the best that can be produced by applying LSD. The 
amplitude variation in the mean line profile due to small 
scale spots is too low to be confidently detected. However, 
Fig. [8] also seems to suggest the presence of smaller ampli- 
tude features at lower latitudes (labelled A, B and C). These 
could possibly be more star-spots at the limit of detectabil- 
ity (they seem to be hinted at in the red tomograms too) or 
they could simply be artefacts of unknown cause. The other 
main feature of the map is the dark region around the Li 
point, best seen in the image at phase 0.5 in Fig. [S] While 
this could be due to a large spot, given the proximity of 
this region to a hot irradiating source (particularly promi- 
nent if the white dwarf is especially hot) it is more likely to 
be due to irradiation. Such irradiation patterns have been 
observed in both high and low resolution ma ps. The irradi- 
ation in RU Peg was previously mapped by iDavevI (119941 ) 
using sodium flux deficit measurements. The map produced 
showed an anomalous region, centred at phase 0.3 on the 
trailing hemisphere and extending some way towards the 
back of the star. Indeed , this feature extended well beyond 
the terminator such that lDavevI (|l994l ') concluded that it was 
highly unlikely to be caused entirely by heating, and instead 
speculated that it might be caused by magnetic activity on 
the secondary star. Our maps show no sign of this feature, 
with th e accretion pa ttern being symmetric. However, the 
maps of iDavevI (119941 ) have no latitude discrimination, so it 
is possible that what was being detected was a superposi- 
tion of symmetric heating and a high latitude spot that was 
located on the trailing hemisphere. 

The right hand panel of Fig. [8] shows the Roche to- 
mogram created using the red-arm data. Again, prominent 



are the off-axis polar spot and irradiation pattern, confirm- 
ing their reality, although in this map the two features are 
rather more blended together. Curiously, there does seem to 
be some non-uniformity in the irradiation pattern with less 
ionisation on the leading hemisphere side of the Li point. It 
is difficult to explain this feature given the lower signal to 
noise of the red profiles and is dismissed as most likely to be 
an artefact. 



Simulating the effects of poor relative flux calibration 

Roche tomography maps relative line fiuxes, and as such, 
the input data are assumed to have an accurate relative fiux 
calibration. In practice, this is achieved by either obtaining 
simultaneous photometric observations, or by observing a 
comparison star on the same slit as the main target. In this 
case we used a comparison star on the slit (as described in 
Section I2.1|l to correct for slit-losses. However, the trailed 
spectrum of the LSD profiles shown in Fig. [5] clearly shows 
some imperfections in the relative scahng of the parent spec- 
tra, despite multiple attempts to remove this effect during 
the data reduction stages. 

In order to assess the potential impact this has on the 
resultant Roche maps, a simulation has been carried out 
using a test image featuring a number of star-spots at a 
range of longitudes. A very low noise trailed spectrum was 
then computed, assuming a negligible intrinsic profile and 
an instrumental resolution of 20kms~'^. Reconstruction of 
the synthetic trailed spectrum, fitting to a x^ of 1.0, showed 
that all the main features present in the test image were 
successfully reconstructed. Next, the trailed spectrum was 
given a 10% error in the relative profile line fiux and re- 
constructions were once again carried out to the same value 
of x^- Again, the reconstructions faithfully reproduced the 
star-spot features, but the stellar surface has an additional 
noise pattern superimposed upon it. It is found that this low 
level noise is removed by slightly under-fitting the data, to a 
X^ of 1.5. Therefore, just as incorrect input parameters have 
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Figure 8. Left - Roche tomogram of RU Pegasi derived from blue-arm ISIS data. The intensity scale is adjusted so that dark scales 
represent regions of lower absorption. The map is presented as the observer would see it at an inclination of 43° except for the central 
image, where the observer is directly above the north pole. Right - Roche tomogram of RU Pegasi derived from red-arm ISIS data. The 
grey scale and viewpoint are the same as in the blue-data map. 



a known effect upon the level of structure in a Roche map, 
imperfect relative flux scaling also produces image artefacts 
that can be identifled. These simulations give us more con- 
fidence that both the irradiation pattern (in general, we do 
not believe the asymmetry apparent in the red data) and 
the high-latitude spot are real features. However, the lower 
resolution of the RU Peg data means that no smaller-scale 
features can be reliably identifled (cf. Fig. [9|. 



6 Hq emission line doppler 
tomography 

Doppler tomography is a powerful technique for mapping the 
emission lines in CVs. Since our RU Peg data cover most of 
the orbit during an outburst, and show prominent Balmer 
emission lines, it seems logical to cons truct a Doppl er tomo- 
gram. Here, we have used the code of lSpruitI (| 19981 ) to map 
the distribution of Ha in velocity space. In this reconstruc- 
tion, the Roche lobe geometry has been set by specifying 
the values for the orbital elements and component masses 
derived earlier. 

It is immediately obvious that there is no strong disc 
emission present in the map; this is probably a result of the 
system being in outburst, and the disc becoming optically 
thick, so that the lines switch from emission to absorption. 
The most striking emission feature is the very bright region 
at the velocity of the secondary star. The brightest pixels 
are found at a position exactly where the inner hemisphere 
of the secondary lies (also marked on the image). This is in 
excellent agreement with our Roche tomograms which show 
symmetric irradiation about the Li point. 



7 CONCLUSIONS 

The Roche tomograms of RU Peg show prominent irradi- 
ation of the inner hemisphere of the secondary star. The 
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Figure 9. Comparison of published multi-line Roche tomograms 
with the RU Peg map in this work. The maps of AE Aqr and 
BV Cen both show the large high latitude spot present in the 
RU Peg map while, interestingly, the provisional V426 Oph map 
does not. Note that the generally smoother appearance of the RU 
Peg map is caused by the fact that the data are of lower spectral 
resolution than for the other three maps, and so cannot reveal 
small-scale features; much of the apparent 'noise' in the other 
maps is probably caused by small spots. 



degree of irradiation found is of no surprise because RU Peg 
was caught near its peak magnitude during outburst. Unlike 
the N a-flux deficit mapping of RU Peg carried out bv lDavevI 
l|l994h . we see no asymmetry in the irradiation pattern, nor 
do we see any evidence for irradiation extending beyond the 
terminator towards the back of the star. 

The Roche tomograms also show a near-polar spot cen- 
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tions (see ISmithl Il995l for a review of the early maps and 
their interpretations). 
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Figure 10. Doppler tomogram of RU Peg in the light of Ha emis- 
sion. Shown are the centres of mass of the primary and secondary 
stars (crosses) and the surface of the Roche-lobe shaped star. The 
ballistic mass transfer stream is also shown. 



tred at a latitude of 82° and covering an estimated 4 per 
cent of the stellar surface. While the resolution of the data 
is too low to place any limit on the overall spot filling fac- 
tor across the star, the presence of a large spot comparable 
to the largest spots seen on BV Cen and AE Aqr suggests 
that RU Peg is also highly magnetically active. As with the 
prominent high latitude spots found on AE Aqr and BV Cen, 
the spot on RU Peg also lies towards the trailing hemisphere. 
(Though we note that, at a latitude of 82°, the spot on RU 
Peg is distinctly more 'polar' than the prominent spots on 
both AE Aqr and BV Cen which lie at a latitude of ~65°.) 
Whether this is a common feature of all CV donor stars and 
indicative of a preferred longitude for magnetic flux emer- 
gence caused by either orbital or tidal effects remains to be 
seen. 

Finally, the spot on RU Peg almost 'merges' with the 
irradiated region. Since spotted and irradiated regions both 
appear as apparent absorption-line flux deficits, we specu- 
late that this may be able to explain the asymmetric irradia- 



tion pattern s found by Na-flux defici t mapp i ng of CV donors 
(see iDav ev 'l994 !: iDavey fc SmithI Il992l : iDavev fc Siiiithl 
119961 : I Cat alan ct al.l 119991 ). Instead of just indicating irra- 
diated regions, could these maps also be affected by the 
presence of large spots which are shifted towards one hemi- 
sphere? With the comparatively low-resolution afforded by 
these maps, it is feasible that large spotty regions could be- 
come blended with the irradiated zone. This inability to dis- 
tinguish between large spots and irradiated regions in such 
methods could lead to maps suggestive of an asymmetric 
irradiation pattern which may appear to extend beyond the 
terminator depending on the location of the spot. This ex- 
planation has the advantage of being able, in principle, to 
explain asymmetries of either sign unlike previous explana- 
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